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d ' Abstract 

^3 : 

O^' Analytical solutions to the time-dependent Shrodinger equation in one di- 

^ i mension are developed for time-independent potentials, one consisting of an 

'tT^ \ infinite wall and a repulsive delta function. An exact solution is obtained 

j_j ■ by means of a convolution of time-independent solutions spanning the given 



Hilbert space with appropriately chosen spectral functions. Square-integr ability 
and the boundary conditions are satisfied. The probability for the particle 
to be found inside the potential well is calculated and shown to exhibit non- 
exponential decay decreasing at large times as t~^. The result is generalized 
for all square-integrable solutions to this problem. 



PACS number: 3.65.Db 



1 Introduction 

Time-dependent quantum mechanics problems are usually addressed using time- 
dependent perturbation theory, adiabatic or sudden approximations as well as several 

*Now at the University of Minnesota, Minneapolis, MN 55455, USA. 
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numerical techniques. It is highly desirable, however, to obtain exact analytical so- 
lutions to given problems, especially in cases when the approximate methods may 
be inadequate to describe detailed aspects of the solutions or when numerical treat- 
ment does not explicitly reveal their mathematical properties at very large times. 
Recently, there has been increasing interest in obtaining exact analytical solutions 
to the time-dependent Schrodinger equation since they can be used to study certain 
physical systems such as quantum dots, Bose-Einstein condensates, unstable compos- 
ite particles and many others. Burrows and Cohen pP have developed exact solutions 
for a double-well quasi-harmonic potential model with a time-dependent dipole field. 
Cavalcanti, Giacconi and Soldati [2j have solved the problem of decay from a point- 
like potential well in the presence of a uniform field and have indicated that, due to 
an infinitely large number of resonances, there may be deviations from the naively 
expected exponential time-dependence of the survival probability. 

The equivalence of exponential decay of a given energy eigenstate with Fermi's 
golden rule when the final density of states is energy-independent and with the Breit- 
Wigner resonance curve has been long known in quantum mechanics jHl E] ■ Dulle- 
mond [Sj has verified this fact in a simple but exactly solvable model and found that if 
final-state energy- dependence is introduced into this model a non-exponential decay 
pattern will dominate at large times. Petridis et al., [Hj have studied a variety of 
systems in which the initial wave function is mostly or entirely set in a finite poten- 
tial well and have observed rich behaviour, including non-exponential decay into the 
continuum. Non-exponentiality for monotonically decreasing survival probabilities 
at large times, though, can be strictly proven only if exact analytical solutions are 
obtained. Specific systems that may exhibit non-exponential decay include systems 
with non-local interactions j7], certain closed many-body systems jHI, quasi-particles 
in quantum dots ^, polarons and non-extensive systems P^!]. 

In this article a method for developing analytical solutions to the time-dependent 
Schrodinger equation is presented. The method is applied to a time-independent po- 
tential, consisting of an infinite wall and a repulsive delta function. In addition to the 
mathematical interest that this potential exhibits it is also applicable to a variety of 
quantum systems undergoing decay. A wavefunction that is an exact solution subject 
to the boundary conditions is obtained and used to analytically calculate the proba- 
bility for finding the particle inside the potential well at any time. A generalization 
of the asymptotic time behaviour to all square-integrable wavefunctions is obtained. 

2 The proposed method 

The method consists of the following steps: (a) The time-independent solutions to 
Schrodinger equation are found subject to the boundary conditions of the problem. 
These are stationary solutions, that is energy eigenfunctions, that span the Hilbert 
space of the given Hamiltonian. (b) Since any finite or infinite, discrete or continu- 
ous linear combination of the stationary solutions (base functions) is also a solution 
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belonging to the given Hilbert space, as long as it is square-integrable, exact analyt- 
ical solutions can be developed by a convolution of the eigenfunctions with energy- 
dependent spectral functions multiplied by the standard oscillatory time-dependence 
of the stationary states. It is, obviously, necessary that the convolution integral over 
the energy converge. This convergence as well as the square-integrability (normaliz- 
ability) of the resulting wave function are verified, (c) The survival probability, i.e., 
the probability for finding the particle inside the potential well is calculated and its 
properties are studied analytically. 



3 Infinite wall and delta-function potential 

The problem to be considered is defined by the one-dimensional repulsive potential, 

V{x)^oo, -oo <x <0 and V{x) ^VoS{x - L), < x < oo, (1) 

with L > and Vq > 0. The steps outlined in the previous section are followed, 
(a) The solutions to the time-independent Schrodinger equation. 



2 dx^ 

(with particle mass m — 1, h — 1 and £^ > 0) are. 



V{x)^e{x)^E^e{x), 



¥^\x) 
¥i\x) 
¥^'\x) 



0, -oo < a; < (region "0"), 

Cism{px), 0<x<L (region"!"), 

C2 sin {px) + C3 cos ipx), L < X < 00 (region "11"), 



(2) 



(3) 
(4) 
(5) 



where p — \/1E and Ci,2,3 are constants in x. These functions obey the boundary 
conditions 



in). 



dx ^ ' dx 



E 



{L)^2Vo^'i\L), 



(6) 
(7) 



while the boundary conditions at a; = are automatically satisfied. The solution is 
not required to vanish at infinity since functions that do not vanish at large x can still 
be solutions to the time-dependent problem. Selecting Ci as the overall normalization 
constant the boundary conditions at x — L yield 



C2 
C3 



Ci 



(2Vo\ 

1 -|- ( 1 sin {pL) cos (pL) 



-cJ^] sin^(pL), 



(8) 
(9) 
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rendering C2 and C3 functions of the energy. The choice of C2 or C3 as the normal- 
ization constant would introduce an energy-dependence in Ci and would effectively 
amount to different spectral functions. 

The obtained linearly independent energy eigenfunctions are orthogonal under the 
inner product 

(^1, ^2) = / i'lix) i>2{x) dx + lim / e-^("-^) iPl{x) 2P2{x) dx, (10) 



e^OJL 



with all wavefunctions in the defined Hilbert space identically vanishing for 2; < 0. 
The orthogonahty relation is 

{^E,'^E')=HE)S{p-p), (11) 



where p — V2E, p — \J2E' and 

w{E) = \\,C2m^\Czm\ (12) 

= l<^i r A V + - 2^)' cos (2pL) + 2pVb sin (2pL)l . 
Ip 

The Dirac 5-function representation used is 

5(x)=lim-— (13) 
(b) The solution to the time-dependent Schrodinger equation, 

can be written as the energy-convolution integral, 

ip(x,t)^ (l)(E)^E(x)e-'^*dE, (15) 
Jo 

with (l){E) a spectral function such that this integral is convergent for all x and all 
t and the resulting wavefunction is square-integrable. Square-integrabiUty of ilj{x,t) 
also requires E to be real. The overall normalization constant is, then, calculated 
from 

/ il^^x.mix.t^dx^l. (16) 
The choice of spectral function 

4>{E) = e-^'^, (17) 

with K a positive constant, produces a convergent energy-convolution integral and a 
wavefunction that is square-integrable even without the presence of the convergence 
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factor that appears in eq. ()10|) . These integrals can be evaluated analytically and in 
closed form. The time-dependent solution is, then, 



0, 



Ci J- X e (K^ + it) 



-3/2 



(18) 
(19) 



iP^''\x,t) 



K-'+^t [K'^ + it) 



-3/2 



eW^)[K^ + it)VQ + 



-it) 



Vo - UVq + x) 



where 



TT 



3/2 



87^3 



+ 



2^3 



TT 



3/2^^2 



2K 



-1/2 



(20) 



(21) 



is the overall normalization factor obtained by means of eq. (jl6|) . 

(c) The probability density p = ip*{x,t)ip[x^t) can be calculated for the interior 
(region I) and the exterior (region II) of the potential well. It is presented in Fig. 1 
at six times starting from t = 0, in increasing order. The initial wave function is not 
entirely localized inside the well. As time progresses the wavefunction spreads and 
tunnels through the potential barrier in both directions. The interference of the wave 
that propagates outwards through the barrier and the wave that is outside creates 
the observed ripples. Inside the well there are no ripples because the wavefunction is 
forced to be odd in x, having a node at a; = 0. The centroid of the probability density 
in region II at t = is always located at 2L, regardless the value of K. 

The survival probability is, then, calculated as 



Pinit) = / ^lj*{x,t)^p{x,t) dx. 

Jo 



(22) 



This yields the closed-form result 



Pinif) — C*! 



TT 



3/2 



87^3 



Erf 



KL 



27rKL 



+ t2 



^K3VK4 + t2 



(23) 



A plot of the survival probability versus time is given in Fig. 2. Pj„(0) is controlled by 
K and approaches an upper limit as L increases at fixed K. For example, this limit 
is equal to 0.9615 for K = 0.1, 0.5 at K = 0.5 and 0.1468 at K = 1.2. It decreases as 
K increases, i.e., as the momentum spectrum becomes sharper. On the other hand 
the decay becomes slower as K increases. The expansion of Pj„ in inverse powers of 
time includes only odd terms with alternating sings. At large times the leading term, 
that has a positive sign, is proportional to t~^, a clearly non-exponential behaviour. 
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Figure 1: The probability density for a potential consisting of an infinite wall and 
a repulsive delta function and using a spectral function that is exponential in the 
energy at six increasing times (from the upper left to the lower right panel, t — 
0.0, 0.3, 0.6, 0.9, 1.2, 1.5). In this plot L = 3, = 1 and X = 1/2. 
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Figure 2: The survival probability for a potential consisting of an infinite wall and 
a repulsive delta function and using a spectral function that is exponential in the 
energy versus time (solid line). In this plot L = 3, Vq = 1 and K = 1/2. The dashed 
line represents the exponentially decaying function, /(t) = aexp(— fitted to data 
points, calculated from the actual solution, in the range t = 2 to 4. The per degree 
of freedom is of order 10~^. 
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4 Corrections to the exponential decay law 

The law governing the decay of physical systems is typically assumed to be a simple 
exponential time-dependence of the number N{t) of the systems that have not decayed 
until time t, i.e, N{t) = A^(0) exp(— At), where A is the decay constant. As mentioned 
earlier this simple law is consistent with the Breit-Wigner curve and Fermi's golden 
rule if the final density of states is energy independent. It refers to the survival 
probability of a given initial energy eigenstate. In the system studied here the initial 
state is not an eigenstate of the energy. If a very large number of systems is assumed 
to be initially described by ip{x, 0) and a system is said to have decayed if the particle 
has exited the potential well, then the number of surviving systems is proportional 
to the probability Pin, i.e., 

N{t) P^t) 



N{0) P,„(0)' 

The differential decay law is 



(24) 



dN = -\{t)N{t) dt, (25) 
where, A is, in general, dependent on time. Substitution from eq. gives 
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Figure 3: The decay parametre A for a potential consisting of an infinite wall and 
a repulsive delta function and using a spectral function that is exponential in the 
energy versus time. In this plot L = 3 and K = 1/2. There is no dependence on Vq. 



In the case studied, eq. yields 

m 



{K^ + t2) [-2ze-'" + Erf (2)] ' 



(27) 



where z = KL/^/K^ + t^. This function is plotted versus time in Fig. 3. 

The decay parameter A peaks later in time and has a smaller maximum value as 
K OT L increase but does not depend on Vq. The expansion of A in inverse powers 
of time includes only odd terms with alternating sings. At large times the leading 
term, that has a positive sign, is proportional to t~^, affirming the non-exponential 
behaviour. At very large times the change of A with time is rather slow. A fit to 
Pin at large times with an exponential curve in a finite time interval (as it is done 
in experiments) gives a very small value of per degree of freedom (of order 10~^) 
so that the distinction between Pj„ at large times and a simple exponential decay 
function is numerically minute (Fig. 2). 



5 Generalization 



Exact, closed-form, analytical solutions to the time- dependent Schrodinger equation 
for the potential consisting of an infinite wall and a repulsive delta function can be 
and have been obtained for other spectral function choices. For example, the choice 



0(E) 



1 



cos 



2E\hL 



{2i 
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yields a square-integrable wavefunction. In the absence of the delta function at x = L 
this would produce an effectively square density pulse at t = located between 
a; = and x = L/2. Due to the actual boundary conditions at x = L this spectral 
function also produces a cusp centered at x = 2L. The survival probability is readily 
expressible in terms of Fresnel sine and cosine integrals. Its asymptotic large time 
behaviour is still t~^. 

A question that naturally arises at this point is whether the asymptotic time 
behaviour can be generalized to other possible solutions to this problem. There is 
a one-to-one correspondence between spectral functions and square-integrable wave- 
functions. This can be seen upon projecting the wavefunction at t = on an energy 
eigenfunction and employing the orthogonality condition of eq. (fTTj) : 

1 r°° 

<P{E) = -— / ^l{x)^>{x, 0) dx. (29) 
w[h) Jo 

Given an initial wavefunction the corresponding spectral function can, in principle, be 
constructed. Schrodinger's time-dependent equation then produces the wavefunction 
at any later (or earlier) time. 

Convergence of the energy convolution integral in region (II) requires that the 
spectral function be finite at E ^ 0. In addition, in order for \l/(x, t) to be square- 
integrable, 0(£') must vanish at large energies. This requirement can be made precise 
by inserting eq. (fT^ into eq. (fTHj) and applying eq. (fTT|) to obtain 

/ w{E)\(P{E)\'^ dE = 1. (30) 
Jo 

Inspection of the function w{E), given in eq. (12), leads to the conclusion that |0(-E)| 
must vanish at infinity faster than 1/^/E due to a constant term in w{E). 

Assuming that 0(£') satisfies the convergence conditions, its contribution to the 
energy convolution integral giving comes mostly from low energies. Then 

at any x in region (I) the wavefunction can be approximated as 

^(^)(x,t) ^ Ci0(O) / V2Exe-'^'dE. (31) 

Jo 

The upper limit of the integration is chosen as follows: the factor exp {—iEt) oscillates 
more rapidly with the energy as t increases. At very large times these oscillations 
eventually lead to a vanishing contribution to the integral. Therefore, the integral 
can be cut off at a point Emax{t) whose first order term in the expansion in powers 
of 1/t is Umax/t, where i/max is constant in t. At low energies (f>{E) is replaced by its 
(finite) value at E = and the sine function is substituted by its argument at a given 
X. Then, the variable change y = Et yields 

ij^^^ ^ Ci0(O) X T'""" e-'y dy. (32) 
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For small Umax the integral is approximately \/2 [(2/3) yflax ~ (2/5) Vmaxl- The 
wavefunction in region (I) is to the first non-vanishing order 

^(^)(x, t) ^ Ci0(O) X M (33) 

where M is a constant and the survival probability decreases with time as t~^. There- 
fore, in order for the wavefunction to be square-integrable, the spectral function must 
be finite at E ^ and decrease at large E faster than 1 / ^/E and then necessarily 
the survival probability asymptotically decreases as the inverse cube of time. 

This argument can be extended to any finite value of x including region (II) 
since the coefficients C2 and C3 are at most of 0(1) for small E. Therefore, the 
integral of the probability density over any finite range of x is finite (even without 
the convergence factor present in eq. (fTUj) ) and it decreases asymptotically as t~^. 

The constant M in eq. ()33|) can be exactly evaluated if (f){E) decreases at large 
E faster than 1/E. Then if (f){E) is analytic in the fourth quadrant of the complex 
E-p\ane the contour integral of (f){E) sin [x\/2E) exp{—iEt) along a closed path, 
consisting of the positive real axis from R to 0, the negative imaginary axis from 
to —iR and a quarter-circle, F, of radius R, is zero. The integral along F is bounded 
by a constant times l/R'^ with R= \E\ and k > 1 and, consequently, vanishes in the 
limit R ^ 00. Then the integration over the real axis gives the same result as that 
over the imaginary axis. The variable change E = —iy with y real, then, yields 

^(^) (x, t) = -iCi (p{-ty) sin (^x ^~2iy^ e"^* dy. (34) 

For large times only small values of y contribute to the integral. The spectral function 
is substantially different from zero only close to the origin and can be replaced by 
0(0) and be pulled out of the integral while the sine function can be approximated 
by its argument in a finite range of x. The remaining integral is easily evaluated as 
a gamma function and gives 

(x, t) ^ Ci0(O) X e-'""^'^^^ (35) 

confirming the earlier result. 

The survival probability, Pi„, calculated thus far refers to the presence of the 
particle inside the potential well. As it has been shown in the previous section the 
spectral function of eq. (fTTjl produces non-zero probability density outside the well at 
t = for i^' > 0. If the "interior" of the well is defined to extend to x much larger 
than 2L (without moving the delta function from x = L) then at t = the probability 
to find the particle "inside" can be arbitrarily close to unity. Specifically the modified 
survival probability Pj„(4L)(t) can be defined by extending the integral of eq. 
to a; = 4L. This integral is then evaluated analytically and plotted in Fig. 4 as a 
function of time. As predicted and verified by an expansion of Pj„(4L) in inverse 
powers of time, its asymptotic time dependence is t^^. An interesting feature of this 
plot is the presence of a step-wise behaviour which can be attributed to interference 
between waves moving in opposite directions. 
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Figure 4: The modified survival probability for a potential consisting of an infinite 
wall and a repulsive delta function and using a spectral function that is exponential 
in the energy versus time (solid line). In this plot L = 3, Vq = 1 and K = 1/2. The 
step-wise behaviour is due to interference of waves moving in opposite directions. 

6 Conclusions and perspectives 

A method has been proposed to solve the time-dependent Schrodinger equation uti- 
lizing the much more easily obtained time-independent solutions for a given Hamilto- 
nian. It has been applied to the case of a potential consisting of an infinite wall and 
a repulsive delta function. An exact, analytical, normalized solution has been ob- 
tained in closed form. The survival probability, which is also analytically calculated, 
exhibits a non-exponential behaviour. At large times it decays approximately as . 
It was shown that this behaviour pertains to all square-integrable wavefunctions that 
are solutions to this problem. To ensure square-integrability the spectral function 
must be finite at £^ — > and decrease to at large energies faster than 1 / \/E. With 
the appropriate choice of spectral functions which, due to linear independence need 
not be the same for waves propagating in different directions, the method could be 
applied to a variety of potentials. It is also of great interest to develop solutions for 
the time-dependent relativistic Dirac equation since these are more appropriate to 
describe the decay of mesons into hght leptons. 
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